MEAN CURVATURE FLOW OF KILLING GRAPHS 
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Abstract. We study a Neumann problem related to the evolution of graphs under mean 
curvature flow in Riemannian manifolds endowed with a Killing vector field. We prove 
that in a particular case these graphs converge to a bounded minimal graph which contacts 
the cylinder over the domain orthogonally along its boundary. 



1. Introduction 

Let M be a (n + l)-dimensional Riemannian manifold endowed with a Killing vector 
field Y. Suppose that the distribution orthogonal to Y is of constant rank and integrable. 
Given an integral leaf P of that distribution, let 17 C P be a bounded domain with regular 
boundary F = dQ. Let : I x Q — )■ M be the flow generated by Y with initial values in 
M, where I is a maximal interval of definition. In geometric terms, the ambient manifold 
is a warped product M = P Xi/^I where 7 = 

Given T G [0, +00), let n : x [0, T) — )■ I be a smooth function. Fixed this notation, the 
Killing graph of u{-,t), t € [0,T), is the hypersurface C M parametrized by the map 

X{t,x) = 'd{u{x,t),x), x^d. 

Notice that this definition could be slightly more general if we suppose that the coordinates 
of X G change with the parameter t G [0,T). To abolish this possibility is equivalent to 
rule out tangential diffeomorphisms of 

The Killing cylinder K over F is by its turn defined by 

(1) K = {§{s,x):sel,xeT]. 

Let be a unit normal vector field along S^. In what follows, we denote by H the mean 
curvature of with respect to the orientation given by A^. We are then concerned with 
establishing conditions for longtime existence of a prescribed mean curvature flow of the 
form 

(2) ^ = - 

(3) X{Q,-) = ^{uo{-),-), 

for given functions tig : — )• M and : 17 — )• M. In order to define boundary conditions 
for the evolution problem ([2]) we consider a function (j) G C°°{T) such that < < 'i 
for some positive constant 0o- Let u be the inward unit normal vector field along K. We 
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where (•, •) denotes the Riemannian metric in M. 

The main result in this paper may be stated as follows 

Theorem 1. There exists a unique solution u : Q x [0,oo) — )• I to the problem 
Moreover, if cj) = and !K = the graphs converge to a minimal graph which contacts 
the cylinder K orthogonally along its boundary. 

Theorem [T] extends Theorem 1.1 in [3J as well as Theorem 2.4 in [2] and Theorem 2.4 
in [T] in a twofold way. The corresponding theorems in [3] and [2] concern evolution of 
graphs in Euclidean space whereas [Tj deals with the case of graphs in Riemannian product 
spaces of the form P x M. Moreover those earlier results hold only for the case when the 
prescribed mean curvature is ^ = 0. Some related results may be also found in fT] and [8]. 

The paper is organized as follows. Section [2] describes the evolution problem in nonpara- 
metric terms. Height and boundary gradient a priori estimates for ([2])-([4]) are presented 
respectively in sections [3] and [H Interior gradient estimates are obtained in Section O Some 
technical computations needed in the body of the proofs are collected in an appendix. 

In Section [6] we prove the following result about the asymptotic behavior of the mean 
curvature flow ([2]) for general (p and IK. 

Theorem 2. Suppose that there exists a solution v G C°°{Cl) of the elliptic Neumann 
problem 



Then the mean curvature flow ©-(jl]) converges to a graph with prescribed mean curvature 
and prescribed contact angle (p. 

In a forthcoming paper the authors prove an existence result for ([S])-®- We also consider 
there another stationary regimen of the mean curvature flow ([2])- ([3]), namely translating 
solitons. 



Since we will consider the mean curvature flow in nonparametric terms it seems adequate 
to describe all geometric invariants as well as their evolution equations in terms of graphical 
coordinates. 

Let x^, . . . ,x" be local coordinates in P. This system is augmented to be a coordinate 
system in M by setting = s, the flow parameter of Y. The tangent space of at a 
point X{t,x), X £ Q, is spanned by the coordinate vector fields 



(5) 
(6) 



div— -7(Vyy,— )=:K ^n ^ 
{u, N) = on dn. 



2. Fundamental equations 



d 



d 



d 
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In terms of these coordinates the induced metric in is expressed in local components by 

(8) gij = aij + ^UiUj, 

where 7 = and (jy are the local components of the metric in P. 

In order to compute the mean curvature of Sj, we fix N as the vector field 

(9) N=^{^Y-^,Vu), 
where Vu is the gradient of « in P and 



(10) W = V7 + IV 



up. 

The second fundamental form of calculated with respect to this choice of normal vector 
field has local components 

(11) a,, = (V^^_^X.A,Ar), 

where V denotes the covariant derivative in M. We then compute 

d - d 

Hence using the fact that the maps x 1— t- 'i?(s, a;) arc isometrics and that the hypersurfaces 
defined by {'&{s,x) : x G P}, s G I, are totally geodesic one concludes that 

It follows from Killing's equation that 

^ ' W W2^ W2j 2W 7"' 

It turns out that aij could be also expressed by 

(13) «« = ^ - ^7(v.v-, A) _ S,(v,r, ^) - ^(V.r, v„). 

Taking traces with respect to the induced metric one obtains the following expression for 
the mean curvature H of the hypersurface St 

^ ^ V WWJ W \ 27 ' ' 



^17) Vyy - , , 



4 J. H. LIRA, G. A. WANDERLEY 

Alternatively one has 

(1^) = r Wt^Jh^ ^7:^7(Vyy, Vn). 

At this point we recall that 

(16) Vj^y = --— y 

and 

1 V7 
what implies that 

(18) {VyY,vu) = -(Vv«y,y) = ^{vj,vu). 

Using this one easily verifies that (jl4p may be written in divergence form as 

(19) d^^-^{V,,Vu)=nH. 
In fact we have 

= U.i ttU U-'Ui-j -u 7,-. 

W);i W 2W^ ' 

It is worth to point out that (llOh is equivalent to 

(20) div— - {VyY, Vn) = nH. 

We conclude that (l2|) may be written nonparametrically as 

du — 

(21) — = l^div— -W:K- 7( Vy y, Vn) . 

Indeed it holds that 

TT r.r ,9X ,du „ d lad. I Ou 

Using (fT4l) one verifies that (f2T]l is equivalent to 

, N du ( \ 27 + |VnP,V7 „ , 



We conclude that the Neumann problem ©-(U]) has the following nonparametric form 

(23) -* = (-^^-^^)-^.-(^ + ^)7^-^-^^ - ^x[0'^) 

(24) n(-,0) = 'uo(-) in x {0} 
with boundary condition 

(25) {N,u) = (p on dnx [0,T). 
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This boundary value problem describes the evolution of the Killing graph of the function 
u{-,t) by its mean curvature in the direction of the unit normal N with prescribed contact 
angle at the boundary. 

The standard theory for quasilinear parabolic equations [Sj guarantees that the problem 
of solving ©-([I]) is reduced to obtaning a priori height and gradient estimates for solutions 
to ([23])- ([25]). 



3. Height estimates 
From now on, we consider the parabolic linear operator given by 

(1 1 \ ■ 

2^ + 2lF2)V-^-^t|F-^--*' 

where e C°°(l^ x [0,r)). 

Proposition 1. For a solution u e C°^(J] x [0,r*]), T* < T, of ([23]) -([25]), it holds that 



Then it follows that 



max = max |ii((0, •)!. 
nx[o,T*] Q. 



_ max |m| < CT* 

nx[o,T*] 



for a given constant C > which depends on T* . 

Proof. First of all we verify that ut is a solution for a linear parabolic equation. Indeed 
one has 

1 1 



Lut = g'^utrj - + ^) (^7, Vut) - utt 

= {9'^u,.j)t - g-iui-j - + (V7, Vut) - Utt 

However since "f = ^{x) in (j22|) and x is independent of t it follows that 



In the same way we have 

(27) Wt = ^{jt + 2u^uk;t) = ^u'^ut.k. 



We conclude that 



Lut = -g'jui-j - r^{V-f,Vu)u''{ut)k + ^^u''{ut)k- 
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Now using the fact that cr*^ = and 7f = we have 

Lut = —(^^^ - ——Wt)ui.i - ^(V-y,Vu)u''(ut)k + —^u^iut)k 
W \ W WW J W^ ' ^ ^"^ w 

,i „,k 



Hence it follows that 

(28) Lut - ^9'\W, - ^){ut)k + ^^{V^,Vu)u\ut)k - ^^u\ut)k = 0. 
Thus fixed T* G [0,r) let (a;o,to) be a point in il x [0,T*] such that 

iif(xo,to) = _max \ut\. 

nx[0,T*] 

Hence we choose a coordinate system adapted to the boundary F in such a way that = v 
at xq. Then, at the point (a;o,to) we have 

Ui-t = ut-i = 

for \ < i < n what implies that 



where we used (j25|) and (p7|) . On the other hand, ([25]) implies that 

(29) Ut-n = Un;t = -{<pW)t = -Hxo)Wt. 

at {xo,to). We conclude that 

{1 - cf>\xo))Un;t = 0. 

However since | (j) \ < 1, it follows that ut-n = what contradicts the parabolic Hopf Lemma 

From this contradiction we conclude that to = 0. Since T* is arbitrary, the conclusion 
follows. □ 



4. Boundary gradient estimates 

Now we will prove a gradient bound for a solution of by applying a modification 

of the Korevaar's technique [1] which appeared formerly in [2j. 

From now on, we consider a non-negative extension d : il — t- M of the distance function 
distp(-,F) satisfying \'Vd\ < 1 in J^. In the same way, we consider a extension of the 
boundary data (j) to the domain which we denote also by (p. Then we define 

(30) r] = e^^'h 

where 



(31) 



h = l + ad-0(Vd, N), 
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where K and a are positive numbers to be fixed later. 

Proposition 2. For a > sufficiently large independent of K and t, if for some t > 
fixed, riW{-,t) attains a local maximum value at a point xq G d^, then W{xQ,t) < K. 



Proof. Let t > be such that 



maxr/VFft, •) = r}W{t, xq) 



for a point xq G T. Hence we choose a coordinate system adapted to F such that = v 

at .tq and 

(32) u\{xq) > and Ui{xo) = 0, for 2 < i < n — 1. 
We have at xq 

(33) = {r]W)i = 7]iW + r]Wi = e^^'iWKuiil - (f^) - 2W(p(pi + Wi{l- (p^)) 
from what follows that 



(34) 


Wi = - 


On the other hand at 


xo we have 


rin = e-^"(ifn„(l- 


^^) + a — (fxf)., 


= e^"(Kn,(l- 


(j)^) + a- 


= e^«(ifu„(l- 


(f>^) + a — 



2(j)(p 



(1-^2) 



Since {riW)n < at a;o it holds that 

> WKUnil - 4>^) + aW- W4>4>n - ^(pUnWn + (t>Un;n + (1 - (t>'^)Wn 

= WKUnO- - <P^) + aW + Wn + (pUn-n + Un<Pn 
= WKUn{l - Cp^) +aW + Wn + cf)Un;n - Wcfx^n- 

On the other hand 

(35) Wn = + — {uiUi-n + UnUn-n) = ^ - r^(pUiWi - (piUi - (pUn-n 

what implies that 

Therefore since 

ul = \Vu\^ -ul = W^ - ^ - (p^W^ = VF2(1 - 02) - 7 



<,2\ 
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we conclude that 

for a given constant C depending solely on 7 and (p. It follows that W{x(),t) < if a is 
chosen large enough and independent of K and t. □ 

5. Interior gradient estimates 

In this section we deduce a global gradient bound using the techniques in [Ij and [2]. 
However the more general context of warped product gives rise to a long list of additional 
terms which require a careful tracking along the calculations. 

In the sequel, we consider the parabolic linear operator given by 

(36) Lv = cf^v,,, - (i^ + ^)l'v, - vt, 
where v G C°°(J7 x [0,r)). 

Proposition 3. For fixed x* <T there exists K > {) sufficiently large so that if 

r]W{xQ,tQ) = max r/VF 
f7x[0,T*] 

for some (xo,to) G x [0,T*], then W{xo,to) < C, for some constant C. 

Proof. We can assume xq ^Vl and to > 0. At a point (xo,to) where v/W attains maximum 
value we have 

(37) r]iW + r/Wi = 
and 

(38) \^^ + ^[^^-^9'^^^^,)<^. 
We conclude that 

1 . 1 „9 „ h 



-Lr] = KLu + -Lh + K'^g'-^muj + 2Kg 



Now we have 
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However 

(iv, ve) = w{AY'^, v^d) + (V V. vd, ^) - AcL-^. 

Therefore 

Thus the expression for Lh in Appendix aUows us to conclude that 

1 97|VnP 
-17] = K'KW + ' 



+^ ( - (iV, Vd) + ^ (iV, V<A>^ + 70(>ii^^, V^d) + ^0(V^ Vd, — ) - 7'^'^^) 
+ '^\A\'^4)e + n^(l)HW{AY^,V^d) + ^(^7 - i(Vd, V7))0(v4y^,y^) 
+^(AV^d,V^,^) + ^(^,v2d)E<A- ■^0(^V^d,X,V7) 

_ 5C) ((AT, V0)0 - a0 - ^(V7, Vd)0 + (V V. Vd, ^)</.) - ni^i/c/) 
-n^Fd + - (2(A^, V(t)) - a) (V w Vd, + -(V Vd, V(/)) 



On the other hand Lemma [3] yields 

^ [lW - ^g'm.Wj) = \A\^ + nHW^AY^, Y^) - nHW^{^,N) 

Now we use the fact that xq is a critical point to r]W. We have 

e^'^iKuih + hi)W = -e^^'hWi. 

what implies that 

-KW'^hNiN' + WhiN' = -hWiN' 

and then 

-K/i|Vnp + WhiN' = -hWiN\ 
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rn;r/3 



V7 



d 



27 



27^ 



^{Vj,N)W + \Vu\'^W{^,N) -W^{AY^,Y^) 
27 27 



and 



hiN' = a9- {V(j), N)9 + (f>aiN'dj - <p{di-jN'N^ - Kaij)N'N^ 



a9 - (V(^, N)9 - (f)W{AY' , V^d) - (/)(V w Vd, — ) + 0k 

W yy 



We then conclude that 



■i^(V7, AT) - 1V«|=(0, N) + W'{AY^, r^) 



Moreover 



w 



!li = Kut + ^ = WK{nH -'K) + ^ 
rj h h 



= nHKW - KWK - ^{nH - J{) N)9 - a9 - -^{"^1, Vd)) 



Then we have 



^[lW - ^g'^WiWj^ = \Af + KnH^ + ^nH - ^nH{V(f>,N)9 - j^nHW{AY^ ,V^d) 
-\{nH - %){lyct>,N)9 -a9- ^(V7, Vd)) - ^nH{V^Vd, ^) + j^nHj-^ 

IV7P 1 



4, ;r^-^'^« + S<-F-^l'^* 



We conclude that 
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where 

A={l + ^)\A\' + ^-^j<P{AY^, V^d) + {{^1- ^(Vd, V7)) {AY^, Y^) 

and 

Vu a 1 \/u 2 

-Jt(V w Vd, — )0 - n-Hd + - {2{N, V0) - a) (V v. Vd, — ) + -(V v. Vd, VcP) 
w W h n ' w W h w 

-i*{V^ V<i, + i^V^M, Vrf) + ni{V/f,, AT)* - i^V^^. ^, |;^) 

However using some standard inequalites we obtain 



Using that VF^ > 7 and choosing a sufficiently large and depending only on n, 7, ^ and k 
we have 

-2'' V ^'/i T^' 27V'' 



27 

>-(e + 2V7^ + ^ + V7|x.Z2|)^ 



Moreover 



^ ^, a a 1 1 1 i^' 



where C is a constant depending on n, 7, ^, d, k and 5C. 
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Hence we obtain 



a K a CK 1 1 1 



Then 



-K^ ' .J > -C — 7 + TT7TT + T77 + TT77 + ^ + TT77T7 + TT77T + TT77T7 + 777 + 777 



/l /l /l /l 



It follows that 



h h h ' J \ h 

+{K + ^ + 1)CW. 
h 

Now suppose that W{xQ,tQ) > 1. Otherwise we are done. In this case we have W < W'^ 
and absorbing the terms with W into that one with W'^ transforms the inequality above 
into 

K^J - -^C-KC-C-^{a + l)(K + 1)C^W^ < (^K^ + K + I + ^{a + l)^C. 

If do = d{xo) then choosing a > l/{C{do)do — l) for some constant C((io) >l/do we obtain 
(1 + a)//i < C{do) what implies that 

- -^C- KC{do) - C(do)) w^<{K^ + K + C{do))C. 



Then for a > ^ max{l, ^/2Cfj} we have 



K^^- KCido) - C((io)) W^<iK^ + K + C{do))C. 



It follows that for K > C{d,)+^C{d,Y+2 jC(d^ ^^^^ j^2j_ KC{do) - C{do) > and 

(39) < CiK^ + K + Cjdo)) 

^ ' - K^-KC{do)-C{do)■ 

This finishes the proof of the proposition. □ 
Theorem 3. There exists a unique solution u :Q x [0, oo) — t- I to the problem ©-([H). 
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Proof. Propositions [H [2] and [3] yield the following global gradient bound 
(40) Wix,t) < Wixo,to)4^^ ^ C^e^^^'^\ 

for {x,t) € n X [0,T*], where Ci and C2 are positive constants and 

M = max \u — uo\. 

nx[o,T*] 

It results that ()23p is uniformly parabolic and then the standard theory of quasilinear 
parabolic PDEs may be applied for assuring the existence of a unique smooth solution to 
(I23D-(I25|). □ 

6. Asymptotic behavior 

In the particular case when u{x,t) = v{x) + Ct, {x,t) £ Cl x [0,T), the initial value 
problem becomes 



Vv - Vv C 

(41) div--7(v.y,-) = :K + - m n 

(42) {ly, N) = (j){x, v) on dQ 

Conversely, notice that if v{x) is a solution of (|41|) -(|42|) then u = v + Ct is a solution of 
()23p which is translating along the flow lines of Y with speed C. 
Now observe that 

S/v - Vf Vf - Vf - Vf S/v 

div- - jiVyY, -) = div- + 7(V v,y, Y) = div- + 7(Vv^r, Y) = div.,^. 

Therefore it follows from divergence theorem that 

(43) / -^+^ = -/ ^W'''^=[ ^^^-)=l 
yi9([o,s]xQ) ^ Ji?([o,s]xr) ^ yi?([o,s]xr) J)?([o,s]xr) 

Since the integrands do not depend on s we have 



(44) 

from what results that 
(45) C 



r 

Jo w 

Since W < Cie'^^^^'^ and < 1 we conclude that 

(46) C < \^\+^JnmMc^^c,MT_ 

Comparing an arbitrary solution of the mean curvature flow with translating graphs 
yields 
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Proposition 4. Suppose that there exists a solution of (|4ip for a given C . Then given a 
solution u{x, t) of i23\) there exists a constant M such that 

(47) \u{x,t) - Ct\ < M 
for {x,t) e 1) X [0, +oo). 

Proof. Consider a solution v of (|4ip . Then consider the functions vi = v + inf (uq — v) and 

V2 = V + sup(uo — v) which are also solutions of (|4ip . By definition we have vi < uq < V2- 
n 

Hence the parabolic maximum principle implies that 

Vl + Ct< u{-,t) < V2 + Ct, 

for t £ [0,T) from what we obtain (I47p . □ 

Theorem 4. Suppose that there exists a solution of (|4ip for C = 0. Then liuit-^oo ut = 0. 
In particular the mean curvature flow converges to a graph with prescribed mean curvature 
"K and prescribed contact angle (j). 

Proof. Since C = Proposition ([4]) implies that 

(48) j <M j \^\ < M\T\, 
for t G [0, r). Also we have 

Therefore 

W dt\Jn Jan J Jn'^W^^ ^' Jn^^W^^ 



an 



n 

It follows that 

f-T 



I f 1^ = - f W{x,T)+ [ W{x,0)- [ u{x,T)cP 
Jo Jn ^ Jn Jn Jan 

^ L " [ L ^^^-^ " [ L ^<^"- ^ ^ 

from what follows that lim w = 0. Since W is bounded then lim ut = 0. This finishes 
the proof of the theorem. □ 

7. Appendix 

In what follows, // and A denote respectively the second fundamental form and the 
Weingarten map of St. Their components are given by 

(50) a, = IHX.^,X.l-):=iAX.^,X.J!-^ 
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Some lemmata will be needed in the sequel. Their content could be also of independent 
interest for other applications. 

Lemma 1. Denote 6 = {Vd,N). The differentials of the functions 6 and h have compo- 
nents given by 

(51) 9i = -aidj + {di-j - Kaij)W 
and 

(52) hi = {a6{ + ^ai)dj - {^{di-j - nuij) + (t)idj)W 
respectively, where k = {'jVyY, Vd) . 

Proof. We have 

QQ Q 

= -{AX,^,Vd) + (AT, ,„.^W) 



'yd — — Vii — — 



Since P is totally geodesic we have 



Moreover we compute 



and 



/Vu - - ,Vu - d . ,Vn , d - 



where we used the fact that Vd] = and that P is totally geodesic. 
Thus we conclude that 



However 



{AX^—^Sd) = ai{X,-^,Vd) =ai{-^ + ujY, Vd) = ajdj = g^'^a^kdj 



Therefore we write 

(53) e, = -g^'^aikdj + {di-j - Kaij)N'j . 

This finishes the proof of the proposition. □ 
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We denote the components of the tensor X*II in P by 
Notice that the covariant derivatives of X*II and // are related by 

However since X^-^ = ^ + UiY we compute 

d d-d - - d - 

_a_X^-— - X^V_a_-— = V _a_—- + Ui kY + mV _a_Y + UkVvi^r- + UiUkVyY 
* 8x1^ ax* ax* axk ax* ax'' ax* 

-V a Tr--(Vtx,V a TT-)^- 
ax* ax* 

Therefore 

5 a - - - 

X -§-^*Tr-: - X*^-S-ir~f = '^i;kY + UiV _a_Y + UfcV a Y + UiUkVyY. 

* dxk ax* ax'' ox^ axk sx' 

Hence using ([12]), (fTUj) and PT|) we obtain 

- 9 9 /I; 1 V7 



VFaifcF + ^UiUkiiVu, Vj)Y + V7) = Wui^Y + ^n^UfcX^VT. 
Hence it follows that 

a X,— -X,V a — ) = MX, — ,— ^ 

' dx^ ' ''^* a^ * 5x* * a^fc 5x* * 5x-? ' 2j 



1 ; UjUfc ; 



We conclude that 



9x*' c^x-? 7 272 



7 27' 



that is, 

(55) VfcOij = Vfc Ojj + -VKaifca^-u/ + - Wajka^ui + -^aja + -^^aia 
Now we use (j55p for computing the Hessian of the function 0. 
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Lemma 2. The trace of the Hessian of in calculated with respect to the metric in S 
is given by 

g^i^ei.^ = -\A\'^e - 2{V'^d, X*II)^ - n{V^H, V^d) - nHW{AY^, V^d) - Ric(Vd, —] 



w 



Proof. Notice that we may write (j53p as 

(56) Oi = -gHiidj + [di-j - Kaij)N^ . 

Hence we have 



= -9'H9^%idj);k + g'^{di;kj + R]kA - fik(Tij)N^ - g'Hdr,j - Kaij){ai - ^k^)- 
However 

g'^{g^%iidj).k = g^^g'^hi-kdj+g'^g'.lbiidj+g'^gHudj.k 

., 5. 1 1 7™ 7"^ 

= 9^ g' (Vfc an + -WaikaY'um + -Waiko'^Um + UiUkaim—ir + uiUkaim-r-^)dj 
7 7 27^ 27^ 

+g'^g'.lhiidj + g'^gHiidj.k 
Hence using Codazzi's equation we obtain 

g'^ig^^biidj).k = g^\nHi + n-WHa^Um + -Wala'^Um + ^fTT^a/^T^ + uiUkal,—,)dj 

7 7 vV^ 27 27^ 



Using that g^^ui = ^^u^ we conclude that 

g'\g'\idj),k = ng^'Hidj - n-W^Hg^'aTNmdj - -W^ g^' 4aT N^dj 

7 7 

+ ^^^kot^dj + g'''gpiidj + g'''g%idj,k 

However we have 
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and 



from what follows that 



g"'{{f%ldj).k = ng^'Hidj - n-W^HaTN^g^'dj - -W^laTNmg'Uj 

7 7 



2 



Therefore 

= -ng^'Hidj + n^W^Ha^Nrng^'dj + ^W^^aia^iV^ff^''^^,- - 2^ 

+g'Hdi;kj + R\kidl - Wj)N^ - g'Hdi;j - Kaij){ai - iVfc^) 
Now using the fact that g^^uj = p^tt* and therefore g^^ Nj = we obtain 



Therefore 

Moreover notice that 

afN' = g^^amiN' = -g>^-^W{AY^ ,X,^) 



and 



aikN^ = -W{AY^,X,^). 
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Similarly we have 



Replacing this above we obtain 



gif^Oi-k = -n{V^H,V^d) -nHW{AY^,V^d) - W{AY^,AV^d) - ^^{AV^d,X^^) 
+W{AY^,AV^d) +j{AY^,Y'^){Vd,^) - \A\'^e - g^^afdj,k 

j 

+9'^{di,kj + R'-jkidi - Kk(Jij)N^ - g'''{drj - Kaij){ai - Nk^) 
Therefore 

gif^Oi-k = -n{V^H,V^d) -nHW{AY^,V^d) - ^-^^{AV^d, X^^) 



+7(^y^, Y^){Vd, ^) - \A\'e - g^'aUr,k 



+g {di;kj + Rjkidi - i^k(^ij)N^ - g'^idi-j - Kaij){al - Nk^) 
However 

Hence we have 

g^'Oi-k = -n{V^H,V^d)-nHW{AY^,V^d) - ^^(AV^d, X.^^) 
+^(Ay^,y^)(Vd, V7) - \A\H - 2g'^g^'di,jaki + ^di-jN'j^ 
+g'''di.kjN^ - Ric(Vd, ^) - ^{N,Vk) + K{nH - j{AY'',Y^)) - ^^(iV, V7) 
This finishes the proof of the Lemma. □ 
Using Lemma [2] we will obtain an expression for Lh. Notice that 

hi-k = Oldi-k - 4>i(^k - (Pkdi - 4>i\kO - (f>0i-k- 

Moreover it holds that 

2<7^Vi^fc = 25^Vi(^V^d,X, A) _ 2g"'dk-i^iN' + iKg^'aki^iX' 
= 2(^V^d, V^cA) - 2g"'dk■l(|)^N' + 2k:^{V^,N). 
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We conclude that 



VF2 

+n(l){y^H,y^d) +n^HW{AY^,V^d) + ^^(t>{AV^d,X^^) 



' * 27 

~^{AY^,Y^){Vd,V^) + \A\^4>e + 2g'^g^Ui,,aki(^ - -^cpd,,^N^^^ 
-g'^di,kjN^ + Ric(Vd, ^)<t> + ^{N, VK)cf> - K{nH - 7(^F^, Y^))<t> 



Now we compute the derivatives with respect to t. We have 

d_ 

= - ( {nH - Oi) , Vd) + {nH - Oi) {N, Vjv Vd) . 

However 



'* = X.-{N,Vd) = {V^eN,Vd) + {N,V^aVd) 



(iV, VatW) = -^(V7, W) + VwW). 



Hence we have 



9t = -{V^'inH - Ji),Vd) + {nH - :k)(-^(V7, Vd) + V w W)). 
Moreover we have 

d - 

(57) dt = ^ > Vd) = {nH - J{) (iV, Vd) = {nH -J{)0. 

Therefore 

/it = a{nH - J{)9 - {nH - Oi){N, V(t))9 + (t){V^{nH - Oi), Vd) 
-</.(n/f - ?{)(-^(V7, Vd) + VwVd)) 

We also compute 

(V7, Vh) = a{Vd, V7) + <^(AV^d, X,V7) - (V<^, V-f)e - ^di-j-f'N^ + k^{N, V7). 
Now we obtain 

g'''di.k = Ad - (V V. Vd, ^) = ~nHd - (V w Vd, ^) 
' w W w W 

V?/ 

g'^dk;i^iN^ = dk-ict^'^N^ - dk;iN^N^N'^i = -(VwVd, V<^) - (VwVd, — )(iV, V(/.). 
' ' ' w w W 

Moreover we have 



and 



5'^V.i;, =A</.-(VwV<^,^) 

w W 
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and 

Therefore grouping and rearranging these expressions we obtain 

Lh = \A\'^(f)e + n(l)HW{AY'^, V^d) + (^7 - ^{Vd, Vj))(f>{AY^, Y^) 

+2{AV^d, V^4>) + 2{A, V2d)s</' - r^4>{AV^d, X^Vj) 

{nH - ((AT, V(/))^ -a9~ ^(V7, Vd)</) + (V w Vd, - riKHcj) 

V?; 

-naHd+ (2{N,V<j)) - a){VvuVd,—) + 2(VvuVd,V(l)) 

-0(V Vd, + 0(V^:K, Vd) + n{VHa, N)cf> - cl>V^di^, ^, + Ric(Vd, ^)</. 
w 27 W W W W 

+^(iV, V.)<^ - + V7) + + ^)(V<^, V7)^ 

Lemma 3. We have 

LW - ^g'^WiWj = \AfW + nHW^{AY'^,Y^) -nHW^{'^,N) -3-f{AY'^,X^^) 
W 27 27 



IV7P 1 



Wt. 

47 

Proof. Notice that 

Wi = -W^{{Vx^^Y,N) + {Y,V^^^N)) 

= -W' {{V^Y, N) + «i(Vyy, iV) - {Y, AX,^)) 



dx^ ' dx'- 



Therefore 



However 
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Hence it follows that 
Hence we obtain 

^g'^w^w, = 4:^^ + ^(V7,iv>(^,iv) + {ay\—1)w^ 
27^ 7 



Now we compute 



However we have 
and 

g'^N^.j = g'^a^kN^^ = -{6i - Nk){a] - N,p) 



27' 



Moreover we compute 



and 



We also have 



ax* 27 7 
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Now we compute 



+Wg^^'^aimN'^-^ + Wg'^N^'^a.ml'^. 
Hence we have 

g'mN%i,j = WN\nVfH + g'^{R{X, — ,X,—-)N,X,—^))+nHW''{AY^,N^X. 



dx^ ' * dx^ ' * dx^ ' * dx^ ' 

+Wh'HAY^,X.±,)(-W{AY^,X.±,)) 



Therefore 



g'mN^bii.j = nWN^VfH - nHW^{AY^ ,Y'^) - W^{AY^,AY^) - \Vu\'^{AY^ ,X^^). 
Moreover 

g'mjN%i = -W^{AY^,^^) + W^{^,N){AY^,Y^) -W^{AY^,AY^) 

2^ z'y 



and 



Wg'mjjki = -Wg'^ia'j - Nj^)au = -\A\^W - ^{AY^ ,X,Vj). 



We conclude that 



g^^Wi-j = \A\'^W + 2W^{AY'^, AY'^) + [nH - , N))W^ {AY^ , Y^) 

+m^\AY'', ^) + \VuWay'^,x,^) + ^(Ay^,x,V7) 

27 472 472 ^ 7 ' ' ^ ^ 27 ' 

-nHW^{^,N)+-iW{VN^,N)-nWN^VfH. 
27 27 



Now 



27 272 
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Hence 

LW - ^g'^WiWj = \A\^W +{nH+l^,N))W^{AY'^,Y'^) 



However 

■7 Z7 Z7 7 Z7 



and 



:|^,iV)^-i(V7.Ar)^W = (|2,iV)', 



and 



= -37(^F^,X.|^) _ W^{^,N){AY^,Y^). 



Moreover we compute 



1 1 x.|V7|%,, 1 



4^2 7 ' 27 ' ' 47 ' 27 

and 

-nM/iV'vf// = ~nW{V^H,N) = -W{V^'K,N) 
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We conclude that 



LW - ^g'^iWi = \ A\'^W + nHW^{AY^,Y^) -3-f{AY^,X^^) 
W 27 



27 4 472 7 ^ 27 ' ' 



However 



and 



7 ^27' ^ 47 w *■ 



7 27 7 27 27 



4 472 ^27'^ 4 472 "^"^ ^4^7'^ ^27'^ 

_ 3|V7lV_LZ2,^)2^. 



4 472 ■ ■ 4 ' 27 
Hence we obtain 



LW - ^g'miWj = \A\'^W + nHW^{AY^,Y^) - nHW^{^, N) - 3j{AY^,X^^) 
W 1^ 27 



- ll^w- - I<l7. + iV) - Vr(V==5(, iV) 



47 w 

This finishes the proof of the lemma. □ 
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